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We study how the ground state of the two-dimensional Ising spin glass with Gaussian interactions 
in zero magnetic field changes on altering the boundary conditions. The probability that relative 
spin orientations change in a region far from the boundary goes to zero with the (linear) size of the 
system L like L'^, where A = -0.70 ± 0.08. We argue that A is equal to d — df where d {— 2) is 
the dimension of the system and d/ is the fractal dimension of a domain wall induced by changes in 
the boundary conditions. Our value for d/ is consistent with earlier estimates. These results show 
that, at zero temperature, there is only a single pure state (plus the state with all spins fiipped) in 
agreement with the predictions of the droplet model. 
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PACS numbers: 75.50.Lk, 05.70.Jk, 75.40. Mg, 77. 

The nature of the ordering in spin glasses below the 
transition temperature, Tc, remains rather poorly under- 
stood. For the infinite range pnoidel, the replica symme- 
try breaking solution of ParisiETcl is generally believed to 
be correct. An important aspect of this solution is that 
the order parameter is a non-trivial distribution, P{q), 
where q describes the overlap of the spin configuration 
between two copies of the system with identical interac- 
tions. The distribution is non-trivial because very differ- 
ent spin configurations occur with significant statistical 
weight. One loosely says that the system can be in many 
"pure states". Monte Carlo simulations on (nroije real- 
istic) short range models on quite small latticesQa, find 
a non-trivial P{q) with a weight at g = which is inde- 
pendent of system size (for the range of sizes studied) , as 
predicted by the Parisi theory. 

An alternative approach, the '^roplet model" , has 
been proposed by Fisher and HuseQ (see also Refs. 
Thermodynamic states and pure states are defined pre- 
cisely by considering correlation functions of spins in a 
region small compared with the system size and far from 
the boundary, and asking whether they change or not 
upon changing the boundary conditions as the (linear) 
system size L, tends to infinity. Each different set of 
correlation functions corresponds to a different thermo- 
dynamic state. The droplet theory, the Parisi theory, 
and some other scftii|a4jios have been studied in detail by 
Newman and SteinEHij. 

By making some plausible and self-consistent assump- 
tions, the droplet theory predicts that the structure of 
pure states is trijvial in short range spin glasses below 
Tc- In zero fieldlHl, trivial pure state structure means 
that any thermodynamic state is a combination of just 
two distinct pure states, related by flipping all the spins, 
which have the same free energy by symmetry. If one 
looks at the whole system, rather than a relatively small 
region far from the boundary, one might note that part 
of the system is in one pure state and the other part in 
the spin-flipped state, with a domain wall between them. 
Hence,a global quantity like P{q) could have a non-trivial 
formll3 even though the structure of pure states is actu- 
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ally trivialEl. 

To unambiguously distinguish between the droplet and 
Parisi pictures it is therefore better to study correlation 
functicms, such as the overlap distribution, in a finite 
regiorJlEI far from the boundary, since the probability 
that the domain wall goes through this region vanishes as 
L — > oo. More precisely, one should investigate whether 
these correlation functions change when the boundary 
conditions are changed. Tq our knowledge, however, this 
has not been done beforalil. 

Here, we perform such calculations numerically for the 
ground states of the Ising spin glass with Gaussian in- 
teractions in two dimensions. Although there is no spin 
glass order at finite temperature in this system, there is 
(complete) spin glass order in the ground state, so one can 
investigate the qjiestion of the number of pure states at 
zero temperatur£3. Two dimensions has the additional 
advantage that there are_ftfficient algorithms for comput- 
ing exact ground statestj'EZi and so quite large sizes can 
be investigated. We find that the probability for the spin 
configuration in the center to change, when the bound- 
ary conditions are altered, goes to zero like as L 
increases, where A can be related to the fractal dimen- 
sion of a domain wall which is induced by the boundary 
condition change. This result shows that there is only a 
single pure state at T = {i.e. a single ground state), 
plus the state with all spins flipped, in agreement with 
the droplet theory. 

The Hamiltonian is given by 



JijSiSj, 



(1) 



where the sites i lie on the sites of an i x i square lat- 
tice with L < 30, Si ~ ±1, and the Jy are nearest- 
neighbor interactions chosen according to a Gaussian dis- 
tribution with zero mean and standard deviation unity. 
Initially, we impose periodic boundary conditions, de- 
noted by "P" . Since the distribution of the interactions, 
Jij, is continuous, the ground state is unique (apart from 
the equivalent state obtained by flipping all the spins). 
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We determine the energy and spin configuration of the 
ground state for a given set of bonds. Next we im- 
pose anti-periodic conditions ("AP") along one direction, 
which is completely equivalent to changing the sign of 
the interactions along this boundary, and recompute the 
ground state. Finally we change the sign of half the bonds 
at random along this boundary, which we denote by "R" . 
Note that the different boundary conditions correspond 
to different choices of the interactions which occur with 
the same probability. Hence they are statistically equiv- 
alent. 

For the smaller sizes, L < 8, we compute the ground 
state by rapidly quenching from a randomly chosen spin 
configuration, and repeating many times until we are con- 
fident that the ground state energy has been found. For 
the two largest sizes, L = 16 and 30, this is impractical 
so instead we use the Cologne spin glass ground state 
serverlla. We repeat the calculation of the ground state 
for the three copies with different boundary conditions 
for a minimum of 2000 samples for each size. 

Next we discuss how to study the dependence of the 
spin configuration on boundary conditions. We consider 
a central block containing Nb = spins, and ask if 
the correlation functions between two spins, i and j say, 
in the block depend on the boundary conditions, i.e. 
whether {S°'S°')t — {S'^Sj)T is non zero for L ^ oo, 
where a and /3 refer to two distinct boundary conditions, 
P, AP, or R here, 5f refers to a spin in the copy with 
the a boundary condition, and (• ■ •)t denotes a ther- 
mal average. We consider even spin correlation functions 
because our boundary conditions do not distinguish be- 
tween states which differ by flipping all the spins. Since 
the difference can have either sign, it is convenient to 
consider its square, i{S^Sf)T - {S'^S'^)Tf. If we sum 
over all the spins in the block, normalize, and average 
over disorder, it is easy to see that this becomes 

A = ((gl)%fe^,)'-2(,f,)^), (2) 

where 

^ Nb 
i—l 

is the overlap between the block configurations with a 
and (3 boundary conditions, and the brackets (• • •) refer to 
both a thermal average and an average over the disorder. 
Eq. can be written as 

A = e [PZX'i) + pU^) '^Pf.M)\ (4) 
where 

P^p{q)-{Hi-0) (5) 

is the probability distribution for the block overlaps. We 
have written these expressions in a general form, valid for 



T > as well as T = 0. Similar arguments can be made 
for correlations of a larger number of spins, which leads 
to expressions like Eq. (^) but with higher moments of 
the overlap distributions. Hence the crucial quantity is 
the difference in the block spin overlap distributions with 
different boundary conditions which occurs in Eq. (Q), i.e. 

AP^p{q) EE P^Jq) + P^^iq) - 2P^,{q). (6) 

If this difference tends to zero as L — > oo then the droplet 
picture is valid. We emphasize that this test does not 
require the size of the block to also become large. 

Specializing now to T = 0, Paai^) is just the sum of 
two delta functions with equal weight at g = ±1, since the 
ground state is unique (apart from overall spin reversal). 
Hence, at T = 0, it is sufficient to investigate the block 
overlap distribution P^/^iq) with a ^ (3. We calculate 
this for a = P, and (3 — AP and R. 

Now we discuss our results, for which we take Lg = 2. 
First of all. Fig. |^ shows data for the root mean square 
difference in ground state energy, 

AE^p^{{El~Elf)''^ (7) 

with the total ground state energy (not the energy 
per spin), for a = P and (3 = AP and R. One sees 
that /S.Ep^AP goes to zero like as L increases, where 
= rOfi^S ± 0.020. This is in agreement with earlier 
workll3il3. The negative value means that large domains 
cost very little energy and so the order in the ground 
state will spontaneously break up at any finite tempera- 
ture, showing that = 0. 
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FIG. 1. A plot of the root mean square ground state 
energy differences AEp^ap and AEp^p for different sizes up 
to L = 30. 

The results for AEp^n are quite different, however, in- 
creasing with L, roughly as L^/^ for large L, rather than 
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decreasing. This difference is easily understood, since the 
defect (i.e. the region where the energy is locally different 
for the two boundary conditions) can be locally removed 
in the P-AP case, by changing the sign of the spins to 
one side of the boundary. The defect will then be a sin- 
gle domain wall somewhere in the sample not necessarily 
near the boundary. However, this can not be done for the 
P-R case and a part of the defect, with an energy which 
one could guess to be ^(''"i)/^ in d-dimensions, will stay 
close to the boundary, in addition to a domain wall which 
could be arbitrarily far away. 

We show some of our data for the block overlaps in 
Fig. ^. The results for the P-AP and P-R overlaps are 
qualitatively similar to each other, with the weight away 
from the peaks at g = ±1 dropping as L increases. 
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FIG. 2. A plot of the block overlaps. 
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and 

Ppii(q), for L = 4 and 16, with block size Lb ~ 2. Note 
that the allowed values of q are 0, ±0.5 and ±1. The left 
hand column is for the P-AP overlap and the right hand col- 
umn for the P-R overlap. The top row is for L — 4 and the 
bottom row for L = 16. The data is normalized so that the 
area under the histograms is unity. 

We characterize this trend by the weight at g = and 
show the results in Fig. ^. For both anti-periodic and 
random boundary condition changes, the weight at g = 
vanishes hkc L^^, where A = -0.70 ± 0.08. This value 
is easy to understand since ^'^-^^(O) is just the probabil- 
ity that the domain wall bisects the block. If the fractal 
dimension of the domain wall is d/ then, generalizing to 
d-dimensions, the probability that it goes through any 
small region is proportional to L~^'^~'^-f\ This immedi- 
ately gives df — 1.30 ± 0.08, which is consistent with 
other estimates for df-. 1.2dJi0.03 by Bray and Moorecil, 
1.34±0.10by Riese* et al.EZland 1.31±0.1 (for a related 
model) by Gingrag23. 
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FIG. 3. A plot of the block overlaps ai q = 
and Ppfl(O), for different sizes up to L = 30 with block size 
Lb equal to 2. 



Earlier calculations have investigated some effects of 
changing the boundary conditions from periodic to anti- 
periodic, usually just the change— in the ground, state 
energy, though Bray and MooreO and RiegertZl have 
also calculated the fractal dimension of the domain wall. 
They obtain a value less than d (as noted above), which 
implies that Pp ap{^) vanishes for large L, as we find 
explicitly here. However, as noted in our discussion of 
Fig. |l|, anti-periodic boundary conditions are special since 
the defect can be locally removed by flipping the spins 
to one side of it. This is why we also investigate random 
boundary condition changes, for which the defect cannot 
be locally eliminated. An important result of our work is 
that the fractal dimension of the domain wall is the same 
in both cases. 

For each configuration of the bonds in the bulk we 
have only studied a single random change in the bound- 
ary conditions. It would be interesting to get statistics 
on a large number of boundary condition changes to see 
if the probability for the domain wall to go through the 
central block obtained by averaging over boundary condi- 
tions for a single large sample is the same as we find here 
by averaging over samples. It would also be interesting 
to investigate boundary conditions which are optimized 
to minimize the ground state energy, i.e. which are cor- 
related with the bonds in the bulk. 

Recently, we have been able to perform calculations 
similar tOjihose presented here for the three-dimensional 
spin glassE3, which has a finite Tc- There too we find 
evidence for a unique ground state. 

After this work was virtually cpjpplete, we became 
aware of related work by Middletonc3. Whereas we start 
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with periodic boundary conditions Middleton takes free 
boundary conditions which allows him to find a ground 
state (for the models studied) in polynomial time, per- 
mitting the study of very large sizes, up to L = 512. In 
this approach one can only perturb the system far away 
from the central region by making it grow bigger. Hence 
there are three relevant lengths: the block size, which we 
call , the size inside which the bonds are not changed 
(which we will call Lmid), and the overall size L. One 
needs Lmid S> and at least some data which also 
satisfies L ^ L^id- Hence the largest sizes L that Mid- 
dleton studies need to be very large. In our work, we 
only have one inequality to satisfy, L 3> , rather than 
two, so the sizes do not need to be as large. Overall, 
the two approaches are complementary and, in our view, 
have similar validity for the two-dimensional spin glass. 
However we believe that our approach is preferable for 
the three-dimensional spin glass, for which there are no 
polynomial algorithms, since it requires smaller sizes. 
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